Discrete Wigner functions and the phase space representation of quantum computers 
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We show how to represent the state and the evolution of a quantum computer (or any 
system with an iV-dimensional Hilbert space) in phase space. For this purpose we use a 
discrete version of the Wigner function which, for arbitrary TV, is defined in a phase space 
grid of 2N x 2N points. We compute such Wigner function for states which are relevant for 
quantum computation. Finally, we discuss properties of quantum algorithms in phase space 
and present the phase space representation of Grover's quantum search algorithm. 
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Wigner functions provide a simple representation of 
the quantum state of a continuous system in phase space. 
This is important when analyzing issues related with the 
classical limit of quantum mechanics In this letter 
we study a useful generalization of the Wigncr function 
to systems with TV-dimensional Hilbert spaces. We use 
this Wigner function to obtain, for the first time, a phase 
space representation of both the states and the temporal 
evolution of a quantum computer. What are the po- 
tential advantages of a phase space representation of a 
quantum computer? It is clear that quantum algorithms 
can be thought of as quantum maps operating in a finite 
state space and are therefore amenable to this represen- 
tation. Whether it will be useful or not will depend on 
properties of the algorithm. Specifically, algorithms be- 
come interesting in the large N limit (i.e. when operating 
on many qubits). For a quantum map this is the semi- 
classical limit where many new regularities arise in con- 
nection with its classical behaviour in phase space. Thus 
unraveling these regularities, when they exist, becomes 
an important issue which can be naturally accomplished 
in a phase space representation. Moreover, this phase 
space approach may allow one to establish contact be- 
tween the vast literature on quantum maps (dealing with 
their construction, semiclassical properties, etc) and that 
of quantum algorithms providing hints to develop new al- 
gorithms. As a first application of these ideas we show 
how Grover's search algorithm |^ can be represented in 
phase space and interpreted as a simple quantum map. 

For a particle in 1 dimension the Wigner function, 
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has the following properties: (PI) W{q,p) is real val- 
ued, (P2) The inner product between states pi and p2 
is Tr{pip2) = 27rn / dqdpWi{q,p)W2{q,p), (P3) The in- 
tegral of W{q,p) along a line in phase space, defined as 
flig -f a2P — as, is the probability density that a mea- 
surement of the observable aiQ -\- a2P has 03 as a result. 
These properties follow directly from the definition (|^). 



However, their origin can be better understood by notic- 
ing that W{q,p) is the expectation value of an operator 
A{q,p)^ known as "phase space point operator": 



W(q,p)^Tr{pA{q,p)). 
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The operator A{q,p), that parametrically depends on 
(q^p), is a symmetrized product of delta functions: 
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Properties (P1-P3) follow directly from simple properties 
of A{q,p). In fact, reality of W{q,p) is a consequence of 
the hermiticity of A{q,p). The inner product rule follows 
from the fact that A{q,p) are complete in the sense that 
rr(i(9,p)i(p',g')) = Siq-q')Sip-p')/27:h. Finally, 
property (P3) is a consequence of the fact that when in- 
tegrating A(q,p) over a line in phase space one gets a 
projection operator: / dpdq5(aiq + a2P — a3)A{q,p) = 
1 13) (03 1 where [03) is an eigenstate of the operator 
aiQ + a2-P with eigenvalue a^. Notice that although the 
above definitions are given in terms of P and Q, only 
unitary exponentials of such operators (i.e, finite phase 
space translations) arc involved in (^). 

To define Wigner functions for discrete systems various 
attempts can be found in the literature. Most notably, 
Woofers [|j proposes a definition that has the properties 
(P1-P3) only when is a prime number. His phase space 
is an A^ X A^ grid (if A^ is prime) and a cartesian product 
of spaces corresponding to prime factors of A^ in the gen- 
eral case. Rivas and Ozorio de Almeida [|j define transla- 
tion and reflection operators relating to the geometry of 
chords and centers on the phase space torus. Bouzouina 
and De Bievre give a more abstract derivation relat- 
ing to geometric quantization. Our approach here follows 
closely Woofers ideas but combine them with the results 
in Q to obtain a Wigner function with all the required 
properties (including the generalization of (P3)). 

Given a basis of Hilbert space (for example, the com- 
putational basis of a quantum computer) = {\n),n = 
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0, .., iV— 1} we can think of it as a discretized coordinate 
basis (with periodic boundary conditions). Its discrete 
Fourier transform Bp — {\k),k — Q,..,N — 1}, related 
to Bx by |fc) = ^^exp(«27rnfc/7V)|n)/\/lV, can be con- 
sidered as the corresponding momentum basis (with the 
same periodicity). For a system with an N dimensional 
Hilbert space it is not possible to define hermitian po- 
sition and momentum operators that generate infinites- 
imal displacements. As Schwinger has shown |^], this 
shortcoming is overcome by introducing unitary oper- 
ators that generate finite displacements. The transla- 
tion operator U generates cyclic shifts in the position 
basis, i.e. U™\n) = \n + m). On the other hand, 
displacement operators in the momentum basis satisfy 
V"^\k) = \k + m) and are diagonal in the B^ basis: 
V"^\n) = exp{i2nmnj N)\n) (analogously, U is diago- 
nal in the Bp basis: f/"|fc) = exp(-i27rfcn/iV)|fc)). The 
commutation relations between U and V take the form 
UV = VUexp{-2in/N] [§. Using this it is easy to 
show that the operator T{m,k) — U'^V'' exp{i7Tmk/N) 
is the discrete equivalent of the symmetrized displace- 
ment exp j-{qP — pQ) in the continuous case. 

With these tools, it is natural to generalize (^ as 

-4(9,^) = ^ E r(m,fc)exp{-z2^fc^}. 

m.k—O 

However, one immediately sees that this definition has a 
serious problem: Thus, for real values of {q,p) the oper- 
ators defined in this way are not hermitian. This prob- 
lem comes as no surprise since the use of an integer grid 
does not allow the half-points needed in (|^). As previous 
works indicate there is an obvious solution to this 
problem which simply consists of defining Wigner func- 
tions in an array of 2A^ x 2N points (lying at half-integer 
values of q,p). This leads to the following definition of 
phase space point operators in the discrete case 

= W ^ nA,A')exp{-^2vriA^-M} 

(4) 

where now q,p are integers in the range {0, 2A^ — 1} 

It is simple to show that these operators have all the 
desired properties and enable us to define the Wigner 
function exactly as in (H) (i.e., the discrete Wigner func- 
tion is the expectation value of the discrete phase space 
point operator First, the operators A{q,p) are Her- 
mitian. Second, they form a complete set. Therefore, one 
can always invert equation (^) and express /5 as a linear 
combination of A{q,p). In fact, the Wigner function are 
the coefficients of such expansion: i.e. 

p = N W{q,p)A{q,p). (5) 
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Using this, it is simple to show that the inner 
product between two states is obtained from the 
corresponding Wigner functions as: Tr{pip2) — 
N^^^Wi{q,p)W2{q,p). It is worth noticing that the 
number of phase space point operators is 4iV^ but that 
there are only iV^ of them that are linearly independent. 
A basis set for such operators is formed by the ones cor- 
responding to points with 0<5,p<7V— 1, belonging 
to the first square subgrid oi N x N points of the origi- 
nal 2N X 2A^ lattice. Thus, the entire Wigner function is 
determined by the values it takes on this first square sub- 
grid. In fact, one can show that W{aqN + q, apN +p) = 
(^—lyqP+'^p'i+Napcrg'^^f^q^p-j^ expression that relates 
the values of the Wigner function on the four different 
N X N subgrids. 

Generahzing property (P3) to the discrete case is less 
trivial. Following Woofers [Q, we can define a line L in 
the phase space grid as the set of points L(ni, 712,713) = 
{(g,p) such that niq + n2P — 713} with < < 2N — 1 
(two lines are parallel if they are parametrized by the 
same integers 7ii and 772 but differ in 773). It is pos- 
sible to show that the above Wigner function satisfies 
the following crucial property: When adding its value 
over all points belonging to a line L(ni,n2,n3) one al- 
ways gets a positive number which is nothing but the 
probability of measuring the eigenstate with eigenvalue 
exp(727r773/iV) of the operator r(77i, 772) (if T(7ii, 772) does 
not have 773 as one of its eigenphases, the sum of W{q,p) 
over L(77i, 772, 773) is equal to zero). This is a consequence 
of the fact that when adding all phase space point oper- 
ators in the line L{ni, 772, 713) one gets a projection oper- 
ator. Thus, 

ii= E i('Z,P) = ^E'^(27.- 773)|7^)(H 

(<?,p)eL 1^ 

where the states are eigenstates of Tn^na with eigen- 
value exp{—i2Trv/N). We have shown that (for even 
values of N) the operator projects onto a subspace 
whose dimensionality is equal to c?^ = Di^/N where 
D]^ is the number of points in the line L for which p 
and q are even (clearly this implies, for example, that 
Dl — if 773 is odd). As an illustration, let us ap- 
ply this to the simplest case: For a line Lq defined as 
q ^ 113, the Wigner function summed over all points 
in Lq is Eiq.,p)eL,Wiq,p) = EpWiq,p) = {q/2\p\q/2) 
if q is even (and zero otherwise). Analogously, consid- 
ering horizontal lines {Lp defined as p — 773) we get 

S(9,P)eL.^(«'^^) = = (p/2|pb/2) if p is 

even (and zero otherwise). These are just two examples a 
general result: this Wigner function always generates the 
correct marginal distributions (this, as in the continuous 
case, is the defining feature of W{q,p)). 

Now, let us discuss the properties of the Wigner func- 
tion of states of a quantum computer. If we consider 
a computational state with p = |77o)(77o| (a position 
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eigenstate), it is simple to show that W{q,p) = SnIq — 
2no)(-l)P('5'-2"'')"/2iV (where zn denotes z modulo N). 
This function, displayed in Figure 1, is nonvanishing only 
along two vertical strips located at q = 2no modulo N. In 
one of these two strips {q ~ 2no) W{q,p) is positive while 
in the other one (g — 2no±N) it oscillates becoming neg- 
ative in odd values of p. These oscillations are typical of 
interference patterns and, in fact, can be interpreted here 
as originating from the interference between the positive 
strip and its mirror image created by the periodic bound- 
ary conditions we are imposing by requiring ciclic behav- 
ior modulo N. The fact that the Wigner function be- 
comes negative in this interference strip is essential to re- 
cover the correct marginal distributions: In fact, adding 
values of W{q,p) along vertical lines gives the probabil- 
ity of measuring a computational state, which is zero for 
all states and 1 for state jng). The Wigner function of 
a momentum eigenstate is entirely analogous to the one 
shown in Figure la but rotated 90 degrees. Another in- 
teresting state is shown in Figure lb where one sees the 
Wigner function of a superposition of two computational 
states. A closed analitic expression is easy to write down 
but the graphic representation is much more eloquent: 
For the quantum state |*) = (|no) + \ni))/y/2, W{q,p) 
is non zero on the two positive strips located at q = 2no 
and q = 2ni and also on a strip located in between these 
two where it oscillates due to interference. The wave- 
length of these fringes depends on the distance between 
the two interfering states and is given by 2N/ (ni — uq)). 
On the other hand W{q,p) is also nonzero and oscillatory 
on the strips obtained from the interference between the 
above three strips and their corresponding mirror images. 
Other states have simple Wigner functions: for the iden- 
tity, W{q,p) is zero everywhere except when p and q are 
both even where it is equal to 1/A'^. 

Temporal evolution of a discrete quantum system has 
also a natural phase space representation. In fact, the 
evolution of p as p' — UpU'^ implies that the Wigner 
function transforms as W'{a) — ZapW{(3), where 
Zap — Tr{AJJApU'^). The unitarity of U implies that 
Za(3 is real and orthogonal (here a and (3 denote points 
in the first N x N phase space subgrid, i.e. a — {qa,Pa) 
with < qa,pa < N). Analyzing quantum evolution in 
phase space is convenient if one wants to reveal proper- 
ties of the semiclassical (large N) regime. The criterion 
for an evolution operator U to have a classical analog is 
simple: This happens if the matrix Zq,^ defines a deter- 
ministic map in phase space. This is so if, for every point 
f3 the matrix element Zap is nonzero only if P — M(a) (in 
this case, the present value of the Wigner function in the 
point a uniquely determines the future value the Wigner 
function at the point /3 = M{a)), This is precisely what 
happens in a classical dynamical flow. In such case we 
can say that the quantum evolution is simply associated 
with the classical map M (i.e., W{a) = W'{M{a))). 



Finding unitary operators with this property is simple 
(but, clearly, generic operators do not satisfy this crite- 
rion). In fact, operators proportional to all phase space 
point operators A{a) are classical in that sense as well 
as all phase space translations Tim, k) and all unitary 
operators associated with linear canonical transforma- 
tions (quantum cat maps Q). Thus, it is simple to show 
that the matrix Za/i associated with U — 2NA{ao) is 
Za(3 — SNia + P—2ao) (corresponding to a classical map 
obtained by a reflection followed by a phase space trans- 
lation). On the other hand, for U = T{m,k) we have 
Zap ~ ^Nia — (3 — 2ao) (which is simply a phase space 
translation by ao = {m,k)). Notice that the factor of 
2 in the above translations is a clear indication of the 
need for introducing a phase space grid of size 2A^ x 27V. 
To the contrary, nonlocal operations are very common: 
Products of Pauli operators acting on different qbits are 
nonlocal in phase space. For example, the matrix cor- 
responding to U = cri"' (where the Pauli operator acts 
on the least significant qbit) is such that Zop = 2/N for 
even values of a being equal to zero otherwise. 




FIG. 1. Left: Wigner function of a computational state of 
5 qubits -black (white) correspond to positive (negative) val- 
ues-. The oscillating strip arises from the periodic boundary 
conditions. Right: Wigner function for a superposition of two 
computational states has two positive strips, an intermediate 
oscillating one and three other strips arising from the inter- 
ference between the above and their mirror images. 

Existing quantum algorithms are a complicated combi- 
nation of classical and quantum operations that, in gen- 
eral, may not have a simple phase space representation. 
However, there are remarkable exceptions: For example, 
the discrete Fourier transform Uft, which plays a major 
role in quantum algorithms has a very simple phase 
space representation (this is not surprising since the very 
notion of conjugate variables relies on Uft)- In fact, 
the Fourier transform is represented in phase space by 
a dynamical map whose action is W'{q,p) = W{—p,q) 
(the Wigner function of the state after applying U ft is 
obtained from the original one by changing the phase 
space point (9,^) into {—p, q), a canonical transformation 
that corresponds to a 7r/2-rotation in phase space). An- 
other important example of a quantum algorithm with 
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a simple phase space representation is Grover's quan- 
tum search whose phase space action is shown in Fig- 
ure 3. In such Figure we display the Wigner function of 
the quantum state of a computer after every iteration of 
Grover's search algorithm. Our system has an = 32- 
dimensional Hilbert space (i.e., the computer has just 5 
qubits) and the algorithm is designed to search for the 
marked item which we chose here to be g = 16 (indicated 
with an arrow in the plot). The initial state is chosen to 
be an equally weighted superposition of all computational 
states which, for the purpose of making the plots more 
visible, we chose to be a non-zero momentum state (the 
usual choice for the initial state in Grover's algorithm is 
a zero momentum state but the algorithm works as well 
with an initial state with nonzero momentum). 



T= 



T= 1 





T= 4 T= 5 

FIG. 2. Phase space representation of Grover's algorithm. 
The computer starts in a pure momentum state and evolves 
(after only five iterations) into a position eigenstate corre- 
sponding to the searched item. 



One clearly observes that the algorithm is very sim- 
ple when seen from a phase space representation. The 
initial state has a Wigner function which is a horizontal 
strip (with its oscillatory companion). After each iter- 
ation VF(q,p) shows a very simple Fourier-like pattern 
and becomes a pure coordinate state at the end of the 
search (in our case, the optimal number of iterations is 
T « 7r\/]V/4 « 5). This representation shows that, as 
a map, Grover's algorithm has a fixed phase space point 
with coordinate equal to the marked item (g = 16 in our 
case) and momentum equal to the one of the initial state. 

We have presented a method to represent both the 
states and the evolution of a finite quantum system in 
phase space. The method has points of contact, and 
also differences, with previous |||,^,||,D and more re- 
cent approaches [0. The main difference that arises 
with respect to the usual continuous infinite dimensional 
case is the fact that finite dimensionality (and periodic 
boundary conditions implicit in modulo N arithmetic) 
produces a characteristic interference pattern between a 
"fundamental" Wigner function and its periodized im- 
ages (clearly illustrated in the Figures). Such fringes, 
that can also be interpreted as arising from the charac- 
teristic aliasing implicit in the discrete Fourier transform, 
are absolutely essential to provide the correct marginal 
distributions from the Wigner function. Finally, re- 
garding potential applications to study quantum algo- 
rithms, we believe that this method provides a novel way 
to analyze their behaviour both in the computational 
and the complementary Fourier transformed basis, thus 
paving the way for a semiclassical analysis of algorithms. 
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